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Abstract. The interplay between fluctuating Cooper pairs and magnetic 
impurities in conventional BCS low-Tc superconductors has been studied through 
measurements of the magnetic field dependence of the fluctuation diamagnetism 
(FD) above Tc in lanthanum with praseodymium impurities. These measurements 
provide a crucial confirmation of our previous observation [Europhys. Lett. 73 587 
(2006)] that in the dilute impurity regime the FD increases almost linearly with 
the concentration of magnetic impurities. This striking effect is attributed to 
a variation due to the presence of the fluctuating Cooper pairs of the coupling 
between magnetic impurities. To describe these results at a phenomenological 
level, we propose a Gaussian Ginzburg-Landau model for the FD which includes 
an indirect contribution proportional to both the impurities concentration and 
the Cooper pairs density. Our approach is able to explain simultaneously the FD 
increase due to magnetic impurities and its decrease with the application of large 
magnetic fields. 
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1. Introduction 

The interplay between magnetism and superconductiv¬ 
ity is still at present an open and interesting subject of 
strongly correlated electron systems [T] . This thematic 
was opened by Matthias et al. [Hi more than fifty 
years ago when they studied the decrease of the super¬ 
conducting transition temperature {T^) of lanthanum 
due to the introduction of magnetic rare earth impuri¬ 
ties, an effect later explained in terms of pair breaking 
by Abrikosov and Gor’kov [1]. In spite of its interest, 
to our knowledge only a few works have addressed the 
interplay between magnetic impurities and the Cooper 
pairs created in the normal state by the thermal agita¬ 
tion (the so-called superconducting fluctuations mm, 
most of them dealing with the fluctuation-induced elec¬ 
trical conductivity (or paraconductivity), Act [7l [H [9] . 

As early suggested on theoretical grounds by 
Schmidt uni, magnetic impurities could enhance 
Act through a change in the relaxation time of 
fluctuating Cooper pairs. Therefore, it could be 
expected that a time independent (static) observable 
like the fluctuation magnetic susceptibility Ay will 
be unaffected by magnetic impurities, apart from 
the effects associated to possible changes in the 
superconducting parameters (for a recent theoretical 
work on this issue see, e.g.. Ref. [H]). This is in 
fact the case of high-Tc cuprates: recent measurements 
|12] have shown that the fluctuation diamagnetism 
of optimally-doped YBa 2 Cu 307_5 is not appreciably 
affected by the presence of magnetic impurities (both 
in the Cu 02 planes or intercalated between them). 
In striking contrast, the measurements in lanthanum 
with diluted magnetic (Pr) impurities presented in 
Ref. [13j have shown an important enhancement of the 
fluctuation diamagnetism, that is observed to be well 
larger than the theoretical results for three-dimensional 
isotropic superconductors laiilllllillllllTlIIHlIIS]. 
The enhancement was found to be proportional to the 
concentration of impurities, being as large as a factor 
of 5 for a Pr concentration of 2 at.%. This effect 
was attributed to a change in the coupling between 
magnetic ions (mediated by conduction electrons) 
due to the presence of the fluctuating Cooper pairs 
[13]. Nevertheless, those results were also affected 
by the experimental uncertainties associated with the 
strong temperature-dependence near (Curie-like) 
of the normal state background, and thus further 
experimental studies of these striking results are highly 


desirable. 

Here we present measurements of the magnetic- 
field-dependence of Ay just above Tc in La and La- 
Pr alloys. The used fields span from 5 x 10“^ T to 
3 T, the latter value being three times larger than 
the upper critical field of the studied samples, where 
fluctuation effects are expected to be negligible IHlEo]. 
The interest of these new measurements is twofold: 
On the one side, the smooth iL-dependence of y in 
the normal state will allow a reliable determination of 
Ay and to check the enhancement with the impurities 
concentration found in Ref. [13]. On the other side, 
these new data will allow to study, for the first time in 
superconductors with magnetic impurities, the finite- 
field (or Prange) regime where fluctuation effects are 
expected to be strongly reduced with respect to the 
zero-field limit [SlEo]. 

2. Experimental details and results 

The samples used in this work (pure La and La-Pr 
alloys) are commercial (Alfa-Aesar and, respectively, 
Goodfellow; 99.9% purity). They were used in the 
previous experiments of Ref. [13], and a thorough 
characterization may be found in Ref. |21j . The 
resulting superconducting parameters are summarized 
in Table 1. The Tc values and the transition widths, 
ATc, were estimated from the temperature dependence 
of the low-field magnetic susceptibility. In the case of 
pure La, the transition temperature is well between the 
ones of pure a-La (~ 5.0 K) and /3-La (^ 6.0 K) 

but no traces were observed of 
diamagnetic steps near these last temperatures. This 
is consistent with a mixing of both crystallographic 
phases at lengths of the order of or smaller than 
the superconducting coherence length amplitude ^(0) 
(~ 20 nm, see Table 1). The ATc values are about 
3% the corresponding Tc values, what will allow us to 
analyze fluctuation effects above Tc down to reduced 
temperatures e = ln(T/Tc) ^ 0.03. In presence of a 
magnetic field, the Tc(iL) shift to lower temperatures 
further increases the accesible temperature range. 
The Ginzburg-Landau (CL) parameter at Tc, k{Tc), 
and the upper critical field (linearly extrapolated to 
T = 0 K), iLc2(0), were obtained from magnetization 
measurements in the mixed state. In view of the k(Tc) 
values, even pure La is a type II superconductor, in 
agreement with earlier results of Pan and coworkers 
[28] . We are not aware of previous Hc 2 measurements 
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in La-Pr alloys, but the results found in pure La are 
consistent with the ones in the literature. For instance, 
in Refs. [531 US [531 [5S] it is found a thermodynamic 
critical magnetic field iiqHc{0) ^ 80 — 84 mT for a-La 
and 110 — 160 inT for ,5-La. By using the relation 
/roiLc(O) = fioHc2{0)/V2K and the data in Table 1, it 
results /ioiLc(O) Ri 95 mT, well within these last values. 
Finally, electron mean-free-paths at low-temperatures 
(estimated from the residual resistivities) are of the 
order or smaller than the corresponding .^(0), so that 
all the samples may be considered close to the dirty 
limit 0 This excludes the presence of non-local effects 
on the fluctuation diamagnetism and, therefore, allows 
an easier comparison with the GL predictions. 

The magnetization measurements were performed 
with a commercial SQUID magnetometer (Quantum 
Design). In order to use all the available volume in 
the sample space, the samples were cut as cylinders 
of 6 mm in height and 6.5 mm in diameter. In 
Fig. [T] we present the magnetic susceptibility of all the 
studied samples as a function of the magnetic field at 
a selected reduced temperature, e = ln(r/rc) = 0.06. 
Even La with 0% Pr presents a paramagnetic behavior, 
probably due to the presence of a small amount (^ 
0.1%) of magnetic rare-earth impurities. The data in 
Fig.H] are already corrected for a small upturn (< 1% 
of the normal state contribution) appearing typically 
below 0.15 T, see Fig. |2])a). This effect is of unknown 
originl^ but it is unrelated to superconductivity 
because it is present in all isotherms up to well above 
Tc (e = 0.5), where superconducting fluctuations 
are expected to be negligible [29]. Therefore, we 
took advantage of its temperature independence and 
corrected the isotherms in the fluctuation region near 
Tc by subtracting the anomaly observed at e = 0.5. An 
example of such a correction is shown in Fig. |5J In any 
case this low-field anomaly has no consequence in our 
analysis, because it focuses in the high-iL region. 

The fluctuation-induced magnetic susceptibility 
Ax may be obtained from the data in Fig. [T] 
by subtraction of the (field-dependent) normal-state 
contribution. This last was estimated for each 
sample by fitting a quadratic polynomial in the 
interval ^ (0.6 — 2)iLc2(0). The resulting background 
contributions are represented as solid lines in Fig. |T] 
The fitting function was chosen because it agrees with 
the data with a good accuracy (the maximum deviation 
in the fitting region is 0.5%) and extrapolates smoothly 
to the lowest helds. Note that to get a better precision 

I In the case of nominally pure La, this is consistent with the 
presence in the sample of a small amount of rare-earth impurities. 
§ It cannot be attributed to the paramagnetic contribution of 
the magnetic impurities, because it should be ff-independent 
when H ^ 0. We also discarded that it is an artifact associated 
to a possible remnant magnetic field in the superconducting coil 
after setting it in persistent mode. 



Figure 1. //-dependence of the magnetic susceptibility just 
above Tc (at e = 0.06). The experimental uncertainty in 
X is of the order of 0.01%, and the corresponding error bars 
are not appreciable. The solid lines represent the normal- 
state susceptibility, determined by fitting a quadratic polynomial 
above 0 . 6 //c 2 ( 0 ), where fluctuation effects are expected to be 
negligible. These data are corrected for a T-independent upturn 
observed at low fields by subtracting the one observed at £• = 0.5 
(see Fig. O. 

in determining the background contribution in the 
low-field region, we extended the lower bound of the 
htting region down to 0.67Lc2(0), neglecting the small 
fluctuation effects that may be present up to ~ Hc 2 { 0 ) 

m- 

Fluctuation effects may be already appreciated 
in Fig. |T] as a x(iL) reduction with respect to the 
background contribution for fields below ^ 0.5iLc2(0). 
The amplitude of this effect is roughly proportional 
to the Pr concentration, as shown in Fig. [3] As 
already commented in the Introduction, this striking 
effect was already observed in measurements of the 
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Sample 
Pr at. % 

Tc 

(K) 

ATc 

(K) 

fJ‘OHc2{0) 

(T) 

m 

(nm) 

i^{Tc) 

i 

(nm) 

0 

5.85 

0.16 

0.8 

20 

4.4 

26.5 

0.5 

5.69 

0.12 

0.9 

19 

5.6 

15.0 

1 

5.51 

0.25 

0.9 

18 

6.0 

11.0 

2 

5.40 

0.18 

1.0 

18 

6.4 

7.5 


Table 1. Summary of the superconducting parameters of the samples studied. The Tc and ATc values were obtained from the 
temperature dependence of the field-cooled magnetic susceptibility measured with = 0.5 mT. //c2(0) and k were estimated 

from the reversible mixed-state magnetization (obtained as the average of the field-up and field-down measurements). The coherence 
length amplitudes were obtained as ^(0) = (0o/27r/io-f^c2(O))^/^. The mean free paths were estimated from the resistivity just above 
Tc by using a Drude model. For details see Ref. m- 



(T) 


Figure 2. Example (for the sample with 1 at. % Pr) of the 
procedure used to correct the x(^) data for the T-independent 
upturn observed at low fields. The line in (a) is a fit of a 
quadratic polynomial to the x{^) data taken at £ = 0.5 (where 
fluctuation effects are negligible), in the H region above the 
upturn. The resulting upturn is then removed from the data 
taken at £ = 0.06, as it is shown in (b). The experimental 
uncertainty in x is of the order of 0.01%. For details see the 
main text. 


fluctuation diamagnetism against temperature in the 
low-field limit m- The present measurements confirm 
those experimental findings, and will allow to study 
the Ay behavior in presence of large reduced magnetic 
helds (in the so-called Prange fluctuation regime). 

3. Data analysis 



Figure 3. Dependence of the fluctuation magnetic susceptibility 
on the impurity concentration in the low field limit {h < 
0.1) and for e = 0.06. These data come from the x{^)t 
measurements in Fig. ^a). For comparison, data from Ref. m 
obtained from xi'^)H measurements are also included. The 
experimental uncertainties in the Ax values are about 2%, and 
the corresponding error bars are smaller than the data points. 


expression is 
= - 
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The Ay data will be analyzed in terms of a 
GL approach with Gaussian fluctuations of the 
superconducting order parameter (GGL approach) 
m i- The direct contribution to Ay (due to the 
Gooper pairs created by thermal fluctuations above 
Tc) was calculated in Ref. m (see also Refs. [HUSO]) 
by using a continuous Landau-level approximation 
(expected to be useful for reduced fields up to h ~ 0.3) 
[19]. This approach includes a total-energy cutoff 
to take into account the existence of a limit (of 
the order of the Cooper pairs size, ^t=ok) to the 
shrinkage of the superconducting wavefunction at high 
reduced magnetic fields or temperatures. The resulting 


£ + q" 


iP 


£ + h + q^ 


-InT 


£ + h + q"^ 


( 1 ) 


2h J \ 2h J \ 2h 

where T and ip are, respectively, the gamma and 
digamma functions, ks is the Boltzmann constant, (pQ 
the flux quantum, and c Ri 0.55 the cutoff constant 
[29] . This expression includes as a particular case the 
low-field (or Schmidt) regime (h <C e), in which Ay is 
field-independent: by imposing h <C e in Eq. m we 
obtain 

Ay««^ = 


tioksmT 


at an. 


3<(>g 




Tc 


( 2 ) 
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Figure 4. Reduced-field dependence of Ax for e = 0.06 in all 
the studied samples. For an adequate comparison these data are 
normalized by ^(0)T. The solid lines represent a fit of expression 
[6], with A ~ —0.5 as the only free parameter. The dashed lines 
are the empirical formula (0.25 -|- 0.6 a:) h~^ (where x is the Pr 
atomic %) evidencing that the extra Ax contribution due to Pr 
ions maintains its proportionality with x also for h ^ 0.3. 

that in absence of cutoff (c -A- oo) leads to the 
classical Schmidt result Ax = —Tr^okBT^{0)/6(j)Q^/e 
M- Equations m and (HI predict the vanishing of 
Ax at e = c, the reduced temperature at which the GL 
coherence length ^(e) becomes comparable to ‘Ct=ok 

m- 

The h-dependence of Ax for all the studied Pr 
concentrations is shown in Fig. |4l To compare the 
results obtained in the different samples, in view of 
Eq. dU these data are normalized by ^{0)T. As 
expected [a [ 20 ], the experimental Ax for pure La 
is in good agreement with Eq. O- However, the 
Ax amplitude increases with the Pr content (almost 
linearly, in view of Fig. El. Moreover, on increasing the 
impurities concentration the h-dependence becomes 
less pronounced. These experimental evidences 
suggest that the magnetic impurities enhance some 
contribution, additional to having a smoother 

h-dependence. As first suggested in Ref. [T3|, a 
natural candidate for the origin of such a contribution 
is the variation of the coupling between magnetic 
(Pr) impurities due to the condensation of Cooper 
pairs. The argument is as follows: in alloys with 
diluted magnetic impurities, the main vehicle for the 
interactions between these impurities is the electronic 
sea. A small variation of the density of normal 
electrons is expected to vary linearly the coupling 
strength between these impurities {e.g., through the 
variation of the effective Fermi vector in the RKKY 
model of electron-mediated interaction between diluted 
magnetic impurities). As the condensation of Cooper 



Figure 5. (a) Tc dependence on the Pr content. The line is a 
fit of the Abrikosov-Gor’kov expression, (b) Linear increase of 
the magnetic susceptibility in the normal state (at a reduced 
temperature e = 0.06), as measured with iiqH = 1 T. The 
experimental uncertainty in x is of the order of 0.01%, and thus 
the corresponding error bars are not appreciable. 

pairs due to superconducting fluctuations will directly 
decrease the density of normal electrons, it may be 
expected that a contribution to the experimental Ay 
will appear in the form 

Ax“ag imp = ( 3 ) 

where ns(h) is the fluctuation-induced superfluid 
density, A is an unknown constant, and x is the 
magnetic impurity concentration (in atomic percents) 
to which Ayamg *mp jg naturally expected to be 
proportional. Similar variations of the magnetization 
associated with the RKKY interaction due to the 
condensation of Cooper pairs and leading to effects 
proportional to ns were proposed previously, although 
not in the context of the superconducting fluctuations, 
by Anderson and Suhl for La-Gd systems in 1959 
m, and more recently for certain borocarbides in 
Refs. [32l|33]. Another reasoning proposed by Bergeret 
et al. [34j based on magnetic moment screening (again 
not related with superconducting fluctuations) also 
lead to contributions to Ay proportional to xns- 

We calculate ns(e,h) for an isotropic 3D super¬ 
conductor on the same grounds as Eq. m (i .e., using 
the GGL approach above Tc with a total-energy cutoff 
and for h < 0.3). For that, we apply to our 0 case 
the same procedure as used in Ref. [29] for h = 0 and 
substitute the h = 0 fluctuation spectrum e -|- ^^(O)A:^ 
by the h > 0 one, e-l- (2n-|- l)h-|-^^(0)fc^. Here k is the 
fluctuation wave vector, kz its component perpendicu¬ 
lar to the applied magnetic field, and n = 0 , 1 ,... the 
Landau-level index. The GGL fluctuation-averaged su- 
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perfluid density is then: 

ns = (l^'l ) °^^e+(2„ + l)/i + ^2(o)yt2 


( 4 ) 


where indicates the statistical and spatial 

average of the squared modulus of the GL wave 
function. Note the appearance in the ^ 

summation of the hT prefactor, that is due to the 
Boltzmann statistical weight and to the degeneration of 
each Landau level. The summation of Eq. o becomes 
then possible using the same continuous Landau-level 
approach and total-energy cutoff procedure as used 
in Ref. [T^] to obtain Eq. (IT|). The result is then 
(employing the same arbitrary dimensionless units as 
in Ref. P5]l: 


ns = 


fioe'^kBT^{0) 

2%"^ J 
e + h + q'^ 


- ^ 


2h 


dq 




c + h + q"^ 
2h 


The solid lines in Fig. 2] are the best fit of 
h) + Ax ns(e, h) 


( 5 ) 

( 6 ) 


to the Ax(/i) data for e = 0.06. The fitting region 
extends up to /i = 0.3, where the above approaches 
are expected to be applicable. As the data in this 
figure are already normalized by ^(0)r, the only free 
parameter is A. The agreement is excellent, in some 
cases even above h = 0.3. The resulting A value (—0.5) 
is reasonably close to the one obtained in Ref. m from 
low-field Ax measurements against the temperature 
[A Ri —0.7), taking into account the important 
uncertainties associated to the background subtraction 
in this last case (mainly affecting the Ax amplitude). 
The data in the region h > 0.3 can be fitted to the 
empirical formula —Ax/C(0)T = (0.25-I-0.6 a;) 
indicating then that the extra Ax contribution due to 
Pr ions maintains its proportionality with x also for 
h > 0.3. 

It is worth noting that Legvold et al. [35] showed 
that the Tc depression with the Pr concentration in 
La-Pr (which is consistent with the Abrikosov-Gor’kov 
theory, see Fig. EKa)) is close to the one of non¬ 
magnetic La-Y and La-Lu alloys. This led them 
to suggest that Pr in the La-Pr alloy could be in 
a singlet state, an that there is little evidence for 
Pr-Pr interactions in these materials. In relation 
with this, in a recent paper Kogan and Prozorov |36| 
calculated that if the upper critical field increases with 
the impurities concentration, the system may not be 
in the purely magnetic scattering regime, and non¬ 
magnetic scattering should also be taken into account. 
However, our observed i7c2(0) variation with the Pr 
concentration is quite similar to the experimental 


uncertaintyljl By this reason, in our simplihed model 
we have ignored a possible contribution coming from 
non-magnetic scattering. Also, the linear increase of 
the magnetic susceptibility in the normal state with the 
Pr concentration (see Fig. |5](b)) clearly indicates that 
Pr does introduce some moderate additional magnetic 
moment in the alloy. Our present results indicate 
thus that even these relatively weak magnetic moments 
(even if perhaps they do not dominate the Tc{x) shift) 
may profoundly affect the diamagnetic contribution 
coming from superconducting fluctuations above Tc- 

The present results contrast with the behavior 
observed in high-Tc cuprates, where the magnetic 
impurities (within the Cu 02 layers or intercalated 
between them) lead to inappreciable effects in the 
fluctuation-induced magnetic susceptibility once the 
change in the Tc value is taken into account. This 
suggests that in these materials the pairs created above 
Tc by thermal fluctuations lead to a negligible effect on 
the interaction between magnetic impurities. 

4. Conclusions 

We have presented detailed measurements of the 
diamagnetism induced by superconducting fluctuations 
in a BGS superconductor (La) under the simultaneous 
presence of magnetic impurities (Pr) and external 
magnetic fields. It is found that magnetic fields 
approaching the upper critical field at T —>■ 0 K 
are very effective reducing the fluctuation effects. 
However, counter intuitively at first, the magnetic 
impurities lead to an important enhancement of the 
amplitude of the fluctuation magnetic susceptibility. 
This enhancement is roughly proportional to the 
concentration of magnetic impurities, and is as large 
as one order of magnitude for concentrations of only 2 
at. % of Pr in La. 

The experimental data are interpreted in terms of 
a change of the coupling between magnetic impurities 
(mediated by conduction electrons), induced by the 
presence of fluctuating Cooper pairs. According to 
this, to analyze the experimental data we proposed 
a Gaussian Ginzburg-Landau model for Ax which 
includes an indirect contribution proportional to the 
impurities concentration, x, and to the Cooper pairs 
density, ns- This approach includes a total-energy 
cutoff to account for the short-wavelength modes 
excited at high reduced magnetic fields, and is 
expected to be applicable up to ~ 0.3i7c2(0). The 
agreement with the experimental data is excellent in 

II See Ref. |21| . where the Hc2{0) values were obtained from 
measurements of the reversible magnetization in the mixed state. 
Due to the highly irreversible nature of the latter in La-Pr alloys, 
Hc 2 { 0 ) was obtained as the average of the field-up and field-down 
contributions. The uncertainty associated with this procedure is 
comparable with the seeming variation of Hc 2 { 0 ) with x. 
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a wide reduced magnetic field interval and for all the 
studied impurities concentration, in agreement with 
previous findings in the low-field limit |13j . It would 
be interesting to further investigate whether fully 
microscopic approaches may account for the amplitude 
of the impurity-induced enhancement of the fluctuation 
diamagnetism. 
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